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I give a quick summary of my proposal for simulating an improvement on quenched QCD with dynamical 
fermions which interact with the gluon configuration only via the topological index of the latter. It amounts to 
include only the topological part of the functional determinant into the measure, thereby absorbing a correction 
factor into the observable. I discuss the prospects of this concept from a study in the massive Af/-flavour Schwinger 
model, where the correction factor is indeed found to be of order 0(1). 



Simulations of lattice QCD have, to date, used 
the quenched and the partially quenched approx- 
imations, where internal fermion loops are omit- 
ted or heavily suppressed by giving the dynamical 
quarks a mass heavier than that of the propagat- 
ing current quarks. Such an approximation turns 
out to be necessary, since presently available up- 
date algorithms with dynamical fermions tend to 
slow down rather dramatically if the fermion mass 
is taken light. From a physical point of view this 
critical slowing down does not come as a sur- 
prise, since it reflects the very nonlocal nature of 
the fermion determinant. The problem with the 
quenched and, to a milder extent, the partially 
quenched approximations, however, is that they 
introduce unphysical degrees of freedom which 
persist in low energy observables and render the 
whole theory sick, if the chiral limit is taken ||l|,|| . 

An alternative path towards full QCD is "To- 
pologically unquenched QCD" |^ . The aim is to 
simulate full QCD, but in an unconventional way: 
The functional determinant is split into a "topo- 
logical factor" , which is included into the measure 
and a "correction factor" which is absorbed in the 
observable. The "topological factor" is defined 
as the determinant of a suitably chosen "stan- 
dard configuration" which shares with the actual 
configuration nothing but the overall topological 
charge. This means that the importance sam- 
pling pays attention only to the first ("topologi- 
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cal" ) factor in the decomposition 

det(^ + m) ^ det(^3td+ m) d^^C^ + m) 

det(^std+ m) 
In the language of full QCD this amounts to a 
statistical guess of the determinant of a configu- 
ration on the basis of its topological index. The 
last factor (i.e. the ratio of the true determinant 
of the configuration to the guess on the basis of 
which it has gotten accepted in the "topologically 
unquenched" Metropolis test) is calculated only 
when a measurement is done (into which it is in- 
cluded). Hence, the number of determinant ratios 
to be calculated is reduced quite drastically - at 
the price of computing for a configuration, when- 
ever possible (see below), its topological index. 

From this outline it is clear that for the concept 
to work two conditions have to hold true. First, 
there must be a correlation between the logarithm 
of the determinant evaluated on a certain config- 
uration and the topological index of the latter. 
Second, the costs (in terms of CPU time) for de- 
termining the topological index must be smaller 
than the savings on determinant evaluation. 

The idea of "topologically unquenched QCD" 
has been launched with a strong focus on the- 
oretical aspects, merely arguing that a correla- 
tion between the functional determinant and the 
topological index (in the continuum theory) is to 
be expected from the Instanton Liquid Model . 
At this time is seems appropriate to do the next 
step and to put it to the (numerical) test. A toy 
theory where the principal issue of "topological 



unquenching" can be studied without having to 
worry about many technical aspects of the im- 
plementation (which do indeed matter in QCD) 
is the multi- flavour Schwinger model (MSM), i.e. 
QED(2) with Nf > 2 light degenerate fermions. 
The physics of this theory resembles in many as- 
pects QCD(4) slightly above the phase transition 
(see refs. [|| and references cited in the second 
of them). Technically, however, there is a differ- 
ence: Even though the gauge fields in the MSM 
(in the continuum-version on the torus) fall into 
topologically distinct classes, there is no precise 
analog of the instanton, i.e. no topologically non- 
trivial object which is simultaneously a minimum 
of the classical action and localized in space-time. 
As a consequence, the "standard configuration" 
in each topological sector cannot be constructed 
artificially but must be chosen as the "most typ- 
ical" configuration in that sector (i.e. the one ^ 
which is in its classical/effective action closest to 
the sectoral average) — cf. strategy (ii) in 

In order to test the "topological unquenching" 
idea, I have chosen to implement the Schwinger 
model with 2 dynamical staggered fermions both 
full and "topologically unquenched", as well as 
the quenched approximation, always using the 
unimproved gauge action ^dass = P X)(l~<^°^(^n)). 
The three theories are compared to each other us- 
ing each time ~ 1500 independent indexed (see 
below) configurations on a lattice of common size 
y=I4 X 6 with periodic/thermal boundary con- 
ditions. The common coupling is P — = 3.0 
and the fermion mass is to = 0.1. These values are 
chosen such that the "pion" (pseudo-scalar iso- 
triplet) has a mass 2.066 5^/^ to^/^c^ 0.3706 
as to fit into the box. The formula used is the pre- 
diction by the bosonized (strong-coupling) ver- 
sion of the theory, which seems adequate, since 
m/g ~ 0.1732 < 1. Note that, unlike in QCD, 
the power by which M-^ depends on to is specific 
for Nf~ 2. Since the operator ^ + to is repre- 
sented by a small (84 x 84) matrix with 5 entries 
per row or column, I have decided to compute the 
determinants exactly, using the routines CGEFA 

^ In practice, it proves useful to relax this condition some- 
what and to construct the "standard determinant" as the 
geometric average of a few determinants of configurations 
which have an action close to the sectoral average. 




Figure 1. Histograms of Sdass, — log(det(_p + TO)) 
[for 2 continuum flavours, after an overall-sliift] 
and S'tot in the full theory with a pair of staggered 
fermions (P = 3.0,m = 0.1) - each for |i^| = 0, 1,2. 

and CGEDI from the LINPACK package. For the 
topological charge, I have implemented both the 
Liischer geometric definition J^geo = ^^^og{Un) 
and the rescaled naive t'nai = ^nai X) '^i'^(^ci) with 
'«nai = 1/(1 - (5'ciass/l^))- A Configuration is as- 
signed an index and used for measurements only 
if the geometric and the naive definition, after 
rounding to the nearest integer, agree (which, at 
/3 = 3.0, holds true for ~ 93% of them). 

The first thing to look at is whether there is 
a correlation between the effective action or its 
components (S'ciass, — log(det(^ -t- m))) of a con- 
figuration and its topological index. The result 
which merely extends work initiated in the liter- 
ature Q is shown in fig. 1 . Only the center of the 
distribution of — log{det{]p + m)) (normalized per 
plaquette) is substantially shifted upwards when 
I J/ 1 increases. The shift between two neighboring 
sectors is about half of the width of the distribu- 




Figure 2. Static quark- antiquark potential and 
distribution of i^nai in the full, quenched and 
"topologically unquenched" (without the correc- 
tion factor) theories, as well as after the correc- 
tion factor (RHS: distribution of its log) has been 
included. LHS: error bars « 2 size of the symbols. 

tion in each of them. Thus there is a positive cor- 
relation between — log(det(^ -\- m)) and and 
the "topological unquenching" idea might work. 

An elementary yet suitable quantity for investi- 
gating the differences between the theories is the 
static quark-antiquark potential, which is shown 
in fig. 2. The full and the quenched potential look 
familiar. The "topologically unquenched" recipe 
generates a sample in between: the raw potential 
(no correction factor being used) shows insuffi- 
cient yet nontrivial screening at long distances. 
The log of the correction factor turns out to be 
sufhciently close to zero to allow for a reconstruc- 
tion of the potential as found with standard (full) 
importance sampling. In terms of CPU time, 
however, there is a big difference: In the present 
example, the standard approach turned out to be 
more expensive by two orders of magnitude. 

The physical reason why "topologically un- 



Figure 3. Potential in the full theory, if only con- 
figurations with h'^Q or v = ±l are evaluated. 

quenched" importance sampling turns out to be 
useful is elucidated by considering the hypotheti- 
cal potentials in the full theory, if only configura- 
tions in a certain topological sector are taken into 
account. As one can see from fig. 3, these "sec- 
toral potentials" agree with the correct potential 
at short distances, but they deviate at larger sep- 
arations. Given the fact that the correct potential 
is, loosely speaking, a weighted average of these 
sectoral potentials, it is clear that it is desirable to 
preserve the relative weight of the different topo- 
logical sectors in the full theory — even if impor- 
tance sampling is done w.r.t. to a different mea- 
sure. Needless to say that this is exactly what the 
"topological unquenching" proposal is aiming at. 

Issues which need more extensive coverage (e.g. 
potential sensitivity on accidental misidentifica- 
tion of V, typical correction factors in larger vol- 
umes, physical reason for the effect seen in fig. 3) 
are hoped to be presented elsewhere. 
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